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fluctuations: general flow equations and application to the 
nematic transition in two-dimensional electron systems 
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We derive renormalization group equations which allow us to treat order parameter 
fluctuations near quantum phase transitions in cases where an expansion in powers 
of the order parameter is not possible. As a prototypical application, we analyze the 
nematic transition driven by a d-wave Pomeranchuk instability in a two-dimensional 
electron system. We find that order parameter fluctuations suppress the first or- 
der character of the nematic transition obtained at low temperatures in mean-field 
theory, so that a continuous transition leading to quantum criticality can emerge. 

PACS numbers: 05.10.Cc, 73.43.Nq, 71.27.+a 
Numerous itinerant electron systems undergo zero temperature quantum phase transi 



tions as a function of a tunable control parameter such as pressure or doping 



l|-|3|. Most 



interesting are continuous transitions associated with quantum critical fluctuations, which 



naturally lead to non-Fermi liquid behavior 



^. However, quantum phase transitions are 



frequently first order, such that critical fluctuations are absent |5|. 

In this letter we explore the possibility that a quantum phase transition which is flrst 
order in mean fleld theory, turns out to be actually continuous by virtue of order parameter 
fluctuations. As a prototypical example we consider the nematic transition jol driven by a 
d-wave Pomeranchuk instability 7|, |8| in two dimensions, which is presently being discussed 
in the context of cuprates [9|], Sr3Ru207 [10|, and other correlated electron materials. Mean 
fleld studies of tight binding electrons on a sq uare lattice indicate that this transition is 

nil] 



typically discontinuous at zero temperature 



13| . A continuous transition occurs only at 



2 



sufficiently high temperature, above a tricritical point. This feature is quite generic since 
the quartic couphng of the Landau free energy describing the transition at zero temperature 
is essentially determined by the curvature of the electronic density of states at the Fermi 
level, but with opposite sign. Due to the logarithmic divergence at the van Hove energy, the 
density of states typically has positive curvature in a broad window around the van Hove 
point, leading to a negative quartic coupling for the order parameter, and hence to a first 
order transition. A continuous transition at zero temperature is found only for judicious 
choices of hopping amplitudes and interactions [l3j . 

A recent renormalization group (RG) analysis of a quantum Landau- Ginzburg- Wilson 
model with 0^ and cf)^ interactions indicated that order parameter fluctuations can alter 
the order of the transition from ffist (in mean field theory) to second order |m]- However, 
that model is well defined only for a positive 0^-interaction, while, for example, the Landau 
energy for the nematic transition discussed above suffers from negative interactions not only 
at quartic order, but rather at any finite order 0^'. The minimum of the Landau energy 
associated with the first order transition and the ultimate increase at large values of cannot 
be accessed by an expansion around = 0. 

In such situations one would like to proceed without resorting to the usual expansion 
in powers of the order parameter field. The functional RG framework offers the possibility 
to approximate the order parameter interaction by a local potential without expanding it 
m powers of HQ. In the following we will derive the functional flow equations for a 
general quantum Landau- Ginzburg- Wilson model for phase transitions in itinerant electron 
systems, and solve them for the specific case of a nematic transition in two dimensions. We 
find that fluctuations can indeed transform the first order transition obtained in mean field 
theory into a continuous one. 

Our starting point is an action for the order parameter of the form 

S[ct^] = \ J^<Pq (^oy^ + ^oq') 0-. + W[0] , (1) 

where is a scalar order parameter field and 0^ with q = (q, w„) its momentum represen- 
tation; Un = 27inT with integer n denotes the (bosonic) Matsubara frequencies. For the 
Matsubara sum and momentum integration we use the abbreviation = TJ2uj„I 
The prefactors Aq and Zq are positive numbers. The dynamical exponent z is restricted to 



3 



values z >2. The potential U[(f)] is a functional of which is local in space and time, 

U[(j)]= I ' dT j d''rU{<P{v,T)) , (2) 



with a (so far) arbitrary real- valued function U (0). The momenta and frequencies contribut- 
ing to S[(j)\ are restricted by an ultraviolet cutoff Aq to the region Aq |q'jl"-2 + Zqc^ < Aq. 

An action of the form Eq. ([1]) can be derived from a fermionic model by introducing 
the order parameter as a Hubbard-Stratonovich field and subsequently integrating out the 



fermions 



17 



18| . Usually is truncated at quartic order, such that only a (quadratic) 



mass term and the local 0^ interaction is kept, which is justified by power counting close to 
a continuous phase transition. Our essential generalization is that we do not truncate U{(j)) 
at all, to be able to deal with cases where an expansion in powers of is not possible. 

We integrate out fluctuations successively by computing the flow of the effective action 
r'^[0] from a functional RG flow equation. The effective action r'^[0] is the generating 
functional for one-particle irreducible vertex functions in the presence of an infrared cutoff 
A. The cutoff is implemented by adding a regulator to the inverse propag ator. r^[0] 
interpolates between the bare action iS[0] for large A and the final effective action r[0] in 
the limit A — )■ 0. Its evolution is given by the exact functional flow equation jiol 

where = d^R^, and r2 [0] = 5^r^[0]/50^. In momentum representation, the trace sums 
over momenta and frequencies: tr = j^. 

The exact effective action is a complicated functional of 0. We resort to an approximation 
of the flow based on the following ansatz: 

r^[0] = + ^q^) 0-. + W^[0] , (4) 

with A-dependent parameters A and Z and a local potential of the form Eq. ([2]) with a 

A-dependent function U{(f)). A classical version of this ansatz, where Matsubara frequencies 
and the dynamical term are absent, has been used previously to analyze classical phase 



transitions 



15 



16| . Inserting the ansatz into the exact flow equation (E]), and evaluating 



the resulting equation for uniform fields, one obtains the flow equation for f/(0), 

dAU{<P) = \ j R\q) G\q- 0) , (5) 



4 



where G^{q; 0) is the regularized propagator in presence of a uniform field 



ICO. 



Aj^^+Zq' + U"{<P)+R\q) 



-1 



(6) 



Here U"{(f)) is the second derivative of U{(f)) with respect to 0, such that the flow equation 
([5]) is a partial differential equation. 

The evolution of the parameters A and Z parametrizing the non-local quadratic term in 
the effective action is determined by the flow of (g; 4>), 



[dp 
1 



R\p) [G\p;<P)rG\p + q; 



-[dpm jR\p) [G\p-cl>)f. 



(7) 



The derivatives dp{(f)) and d^U{(f)) correspond to local 3-point and 4-point vertices, respec- 
tively. We extract the Z-factor from the two-point vertex as Z = ^ Aq [Fg (f/; 0o)] |g_o > 
where 0o is the position of the global minimum of f/(0). The Laplacian Aq is defined at 
fixed R^, that is, it does not act on the regulator R^{p + q) in G^{j) + q) on the right hand 
side of Eq. Q. The flow of the A-factor could also be extracted from F^, but it is of minor 
importance |2o| and will be discarded here, that is, we set A = Aq. 



As a regulator we choose a Litim-type [2l| cutoff function 



R\q) = [Z (A^ - q2) - A|qp-1a;„|] 6 [Z (A^ - q^) - A\qMu;^\] , 



(8) 



which has the convenient feature that the regularized progagator becomes momentum- and 
frequency-independent in the integration region, that is, (^^(g; 0) = [ZA"^ + f/"(0)]~^ for 
Alqp"^!^^! -|- Zq^ < ZA^. The momentum integrations on the right hand side of the flow 
equations can be carried out analytically. The frequency summation in Eq. is performed 
numerically. 

The flow of Z receives significant contributions only near a continuous phase transition at 
finite temperature, where non-Gaussian classical fluctuations generate an anomalous scaling 
dimension (given by the logarithmic derivative of Z with respect to A) 20j. Note that 
quantum critical points at T = are Gaussian for d + z > A [17]. We can therefore neglect 
the contributions from 7^ to the flow of Z, and obtain 



dlogZ 
dlogA 



-T[U"'{<j>o)? 



A" 



{27tY Z3(l-f5)5 d 



2(1 + 5) 



d + 2 



(9) 
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where 5 = U" {(f)o) / [Z A'^) and Sd-i = /T{d/2) is the area of the (rf-l)-dimensional unit 
sphere. 

We now apply the flow equations to a specific model system, featuring a nematic transition 
driven by a d-wave Pomeranchuk instability. The model describes tight binding fermions on 
a square lattice with an attractive d-wave forward scattering interaction. Its Hamiltonian is 
given by Q 23 1 

^ = Xl^knk + — ^ /kk'(q)^^k(q)^k'(-q) , (10) 

k k,k',q 

where nk(q) = Xlo- '^k-q/2 o-'^k+q/2 o- ^^"^ ^ number of lattice sites. For nearest and 

next-to-nearest neighbor hopping the kinetic energy reads ek = — 2t(cos kx + cos ky) — 
At' cos kr COS ky. It s van Hove (saddle) points are situated at k = (vr, 0) and (0,7r). The 
interaction has the form 

/kk'(q) = -fi'(q)c?kc^k' , (11) 

where lik = cosk^ — cos ky is a form factor with d^^-y^ symmetry. The coupling function 
5'(q) > has a maximum at q = and is restricted to small momentum transfers by a 
momentum cutoff. For sufficiently large g = g{0) the interaction drives a d-wave Pomer- 
anchuk instability leading to a nematic state with broken orientation symmetry, which can 
be described by the order parameter 

= fE^k(nk). (12) 

k 



The mean-field solution of the model ffTOj) has been analyzed in a series of articles llNl3|. 
In the plane spanned by the chemical potential and temperature a nematic phase is formed 
below a dome-shaped transition line Tc(/i) with a maximal transition temperature near van 
Hove filling. The phase transition is usually first order near the edges of the transition line, 
that is, where Tc is relatively low, and always second order at the roof of the dome. 

Introducing an order parameter field via a Hubbard- Stratonovich transformation, inte- 
grating out the fermions, and keeping only the leading momentum and frequency depen- 
dences for small q and aJ^i/lql leads to an action S[(p] of the form Eq. ([1]), with z = 3 
and 

f/(0) = I - 2T I ^ In (1 + e-^--^^-'^)/^) . (13) 

The coefficients of a Landau expansion of this potential, f/(0) = + |f 0"^ + . . . , are given 
by 02 = g'^ - A^2(/i) and a2i = -|^A^2z(/x) for / > 2, where A^2z(/x) = x Zlk c^k'^(ek - /^) is 



6 



the density of states at the Fermi level with a d-wave form factor [13]. Since A^2/(/^) diverges 
logarithmically at van Hove filling, a2i is typically negative for I > 2. 

The coefficients Aq and Znin S[(j)] are related to the d-wave particle-hole bubble and 
the coupling function g{q) 23|]. Note that Aq and Zq are not involved in the mean- field 
solution of the model Eq. (ITOl) . Aq is associated with the Landau damping term in the 
bubble, and depends actually on the orientation of q. We will neglect this q-dependence 
and insert a constant Aq corresponding to an angular average. While the bubble is fully 
determined by the kinetic energy ek, the choice of Zq remains to a large extent arbitrary 
since it also depends on the momentum dependence of 5'(q). Since momenta q in the action 
S[(f)] correspond to momentum transfers in the underlying fermionic model, the ultraviolet 
cutoff Aq in S[(f)] limits momentum transfers in g{q) to |q| < Ao/\/Zq. 

We now present results for the nematic transition as obtained from the functional RG. 
The following parameters are the same in all quantitative plots: t = 1, t' = —1/6, g = 0.8, 
Aq = 1, and Zq = 10. In Fig. 1 we show an exemplary plot of the evolution of the flowing 
effective potential U{(f)) for A ranging from the ultraviolet cutoff Aq = e^^ ~ 0.37 to the 
final value A = 0. The initial (mean-field) potential has a minimum at (po = 0.112. The 



-0.8932 



-0.8936 



-0.894 



-0.8944 



-0.8948 



K=0.22 
A=0.11 
A=0.05 
A=0.00 




0.05 



015 



0.2 



FIG. 1: (Color online) Effective potential U{(j)) for various values of A between Aq = e ^ and 0, at 
/i = -0.78 and T = 0.05. 



final potential exhibits spontaneous symmetry breaking with an order parameter 0o = 0.102. 
Fluctations shift 0o toward a slightly smaller value compared to the mean-field solutions. 



Note the fiat shape of ?7(0) for < 



of the grandcanonical potential 



24 



at A = 0, which is imposed by the convexity property 



251]. The final value of 0o can be obtained either by 



7 



following the sequence of minima 0q until A — )■ 0, or by determining the point on the 0-axis 
where the fiat shape in the final potential U{(f)) terminates. 

In Fig. 2 we present the transition line between normal and symmetry-broken phases 
for two choices of Aq. Compared to the corresponding mean-field result, the transition 
temperature is suppressed, with a larger reduction for larger Aq (corresponding to a larger 
phase space for fluctuations). While the phase transition is of first order at low temperatures 
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FIG. 2: (Color online) Critical temperatures versus chemical potential for Aq = (larger dome 
with dots and crosses) and Aq = 1 (smaller dome with dots). The mean- field transition line is also 
shown for comparison. 
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2nd order transition 
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2nd order (MFT) 



in mean- field theory, and also in the presence of fluctuations with a cutoff Aq = e~^, for 
Aq = 1 the transition is continuous down to T = 0, leading to quantum critical points /ic 
at the edges of the nematic dome. In Fig. 3 we show the order parameter 0o at T = as 
a function of /i near the left edge of the transition lines in Fig. 2 for Aq = and Aq = 1. 
For the case Aq = the jump of (j)Q indicates a pronounced first order transition. By 
contrast, for Aq = 1 the order parameter rises continuously, and the shape of (t)o{n) exhibits 
the expected behavior of a Gaussian quantum critical point. 

We note that the mean-field phase diagram is fully determined by the parameters t, t'. 



robust 



the 


size 


12. 


13| 



13|. The parameters Aq, Zq, and Aq determine the strength and phase space of 



fluctuations. The suppression of and of the first order lines is stronger for larger Aq and 
smaller Aq and Zq. 

In summary, we have presented a renormalization group method which allows us to treat 
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FIG. 3: (Color online) Order parameter (j)Q as a function of /x at T = near the left edge of the 
nematic dome for Aq = e^^ and Aq = 1. 

order parameter fluctuations near quantum phase transitions in cases where an expansion 
in powers of the order parameter (f) is not possible. We have derived flow equations for the 
full effective potential U (0) in a local approximation. As a first application of the method, 
we have analyzed fluctuation effects on the nematic transition driven by a d-wave Pomer- 
anchuk instability in two-dimensional itinerant electron systems. It turned out that order 
parameter fluctuations suppress the first order character of the nematic transition obtained 
at low temperatures in mean-field theory, such that a continuous transition associated with 
quantum criticality can emerge. Our flow equations may also be applied to a ferromagnetic 
transition with Ising symmetry and to other Pomeranchuk instabilities. 
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